Short-time motion of Brownian particles in a shear flow 
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The short-time motion of Brownian particles in an incompressible Newtonian fluid under shear, 
in which the fluid inertia becomes important, was investigated by direct numerical simulation of 
particulate flows. Three-dimensional simulations were performed, wherein external forces were in- 
troduced to approximately form Couette flows throughout the entire system with periodic boundary 
conditions. In order to examine the validity of the method, the mean square displacement of a sin- 
gle spherical particle in a simple shear flow was calculated, and these results were compared with a 
hydrodynamic analytical solution that includes the effects of the fluid inertia. Finally, the dynam- 
ical behavior of a monodisperse dispersion composed of repulsive spherical particles was examined 
on short time scales, and the shear-induced diffusion coefficients were measured for several volume 
fractions up to 0.50. 

PACS numbers: 82.20. Wt, 82.70.-y, 05.40.-a, 83.50.Ax 
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I. INTRODUCTION 

The short-time motion of small particles fluctuating 
in a Newtonian fluid is strongly affected by fluid inertia, 
that is, the vorticity of the host fluid surrounding the dis- 
persed particles. If a dispersed particle accelerates due to 
Brownian forces, it affects the motion of the host fluid in 
the neighborhood of the particle, while the vorticity gen- 
erated by the particle's motion then affects the motion of 
the same particle. These effects are referred to as mem- 
ory effects and have an important role in the dynamical 
motion of a dispersion on short time scales [H, 0, El . 

The vorticity diffuses away on a kinematic time scale 
t v = a 2 jv, where a is the radius of the spherical par- 
ticle and v is the kinematic viscosity of the host fluid 
(~ 10~ 6 s for 1 /im in water). When the vorticity dif- 
fuses away much faster than the particle's motion, i.e. 
t v <C tb = M /Qitrja where M is the mass of a Brownian 
particle and rj is the shear viscosity of the host fluid, the 
motion of a Brownian particle is well approximated by 
the normal Langevin equation, which is valid for strong 
damping (Reynolds number Re — > 0) or long time scales; 
however, this equation, wherein the effects of the fluid 
inertia are ignored is not applicable to a dispersion com- 
posed of neutrally buoyant particles since tb is compara- 
ble to r„. For a complete understanding of the short-time 
motion of a dispersion, the inertias of the particle and the 
host fluid cannot be neglected. 

One way to account for memory effects is to simulta- 
neously resolve the fluid motion with the particle motion 
as a boundary condition to be satisfied. We refer to this 
approach as the direct numerical simulation(DNS) ap- 
proach. Within the DNS approach, various numerical 
methods have been developed [J, [H, @, 0, H, [1] , and the 
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power-law decay behavior in the velocity autocorrelations 
of a free Brownian particle has been accurately repro- 
duced. This slow relaxation of the correlation behavior 
is one of the main features of memory effects. Although 
most of these methods have been applied to dispersions 
composed of free Brownian particles at thermal equilib- 
rium, dispersions under flows that are far from equilib- 
rium have not been examined in detail, even for the sim- 
ple case of a single Brownian particle in a shear flow on a 
short time scale, r„. Furthermore, most numerical meth- 
ods used for concentrated dispersions under shear, which 
are widely used for measuring rheological properties, are 
limited to a Reynolds number of zero, and the short-time 
motions of the dispersed particles in concentrated disper- 
sions cannot be correctly tracked. 

Recently, we have developed a numerical method, 
known as the "Smoothed Profile Method (SPM)" 0,[H, 
for the DNS of particulate flows. Its computational ac- 
curacy and efficiency have been examined carefully by 
Lio et al. for several flow problems [12]. The SPM has 
been applied to a dispersion com pos ed of free Brownian 
particles at thermal equilibrium [l3l . [l4| . We have also 
succeeded in reproducing the power-law decay behavior 
in the translational and rotational velocity autocorrela- 
tions of a Brownian particle, and these results agree well 
with hydrodynamic analytical solutions for a free Brown- 
ian particle in an infinite fluid that accounts for memory 
effects. 

In order to simulate dispersions in non-equilibrium 
conditions on short time scales in which memory effects 
become significant, we have modified the SPM. The pri- 
mary objective of the present work is to accurately ex- 
amine the short-time motion of Brownian particles in a 
simple shear flow by using the modified SPM. In this pa- 
per, we first present the modified SPM, in which exter- 
nal forces are introduced to impose a shear flow into the 
system. In order to validate the method, we next com- 
pare the numerical results for the mean square displace- 
ment (MSD) with a hydrodynamic analytical solution of 
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a Brownian particle in simple shear flows that account 
for fluid inertia. Furthermore, we apply our method to 
a dispersion composed of many spherical particles under 
shear. The MSD in the vorticity direction is then calcu- 
lated for several volume fractions, and the time evolution 
is discussed. Finally, the shear-induced diffusion coeffi- 
cients are measured from the long-time behavior of the 
MSD, and the dependence of the determined diffusion 
coefficients on the volume fraction is examined. 



II. SIMULATION METHOD 

Let us consider a monodisperse dispersion of repul- 
sive spherical particles in a Newtonian host fluid. The 
dispersion is subjected to shear by an external force. 
The position of the ith particle is Ri, the transla- 
tional velocity is V*, and the rotational velocity is fi^. 
The velocity and pressure fields of the host fluid are 
v(x,i) and p(x,t), respectively. These field quanti- 
ties are defined on three-dimensional Cartesian grids; 
x E [0,L x ][—Ly/2,L y /2][0,L z ]. In order to distinguish 
the particle and fluid domains on the grids, a smoothed 
function <p{x,t), which is equal to I in the particle do- 
mains and in the fluid domains, is introduced. These 
domains are separated by a thin interfacial domain of 
thickness £. The length unit is taken to be the lattice 
spacing A, and the time unit is p/A 2 /?y, where pf de- 
notes the density of the host fluid. The time evolution of 
the ith dispersed particle with mass Mi and moment of 
inertia Zj is governed by Newton's equations of motion: 

M i V i = F? + Ff + GY, Ri = Vi, (1) 
Ii-ni = Nf +Gf, (2) 

where F t H and NP are the hydrodynamic forces and 
torques exerted by the host fluid on the particle [l(J EH • 
Ff is a repulsive force that is employed to prevent par- 
ticle overlaps, and a truncated Lennard- Jones potential, 
V{ra) = 4[( C j/r ii ) 36 - (<r/ry ) 18 + 1/4 for r v] < 2 1 / ls a or 
V(rij) = 0, is adopted in this work. Here r,j = \R% — Rj\, 
and a = 2a represents the diameter of particles. GY and 
are random forces and torques, respectively, due 
to thermal fluctuations. These fluctuations are intro- 
duced as white noise with a zero mean and correla- 
tions (GY(t)GY(0)) = a v I8{t)5 t] and (G?(t)Gf(0)) = 
a n IS(t)5ij , where a v and a n are numerical parameters 
that control the translational and rotational particle tem- 
peratures, namely, T v and T Q [U, El- 

The angular 

brackets denote taking an average over an equilibrium en- 
semble. The temperatures are determined by the follow- 
ing procedure; First a single Brownian particle at ther- 
mal equilibrium is simulated with fixed ay and an . Then 
the translational and rotational long-time diffusion coef- 
ficients are obtained from the simulation. By comparing 
these diffusion coefficients with Stokes-Einstcin diffusion 
coefficients, DY = ksT v /Sirr/a for the translational mo- 
tion and Dg = hsT /8irrja 3 for the rotational motion, 
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FIG. 1: Schematic diagrams of a zigzag velocity flow. The ar- 
rows represent the velocity vectors of the host fluid in the flow 
direction. The x, y, and z axis represent the flow, velocity- 
gradient, and vorticity directions, respectively 

we finally can determine the temperatures. Since both 
ay and an are chosen to satisfy fcsT y = ksT n , the 
temperatures are simply written as fcgT throughout this 
paper. 

The time evolution of the host fluid is governed by the 
Navier-Stokes equations: 

p f {d t v + vVv) = V-<t + Ps Mv + Pff shear , (3) 
er = -pI + T]{Vv + (Vt>) T } (4) 

where the incomprcssibility condition, V • v = 0. 
f shear (x,t) is an external force field that is introduced 
to enforce the following velocity profile, 

(j(-y-L y /2), (-L y /2<y<-L y /4) 
v x (v)=\iy, (-L y /4<y<L y /4) (5) 

(i(-y + L y /2) (Ly/4<y<Ly/2) 

where 7 denotes the shear rate of the imposed flow and 
y denotes the distance in the velocity-gradient direction. 
This velocity profile, schematically depicted in Fig. [I] en- 
ables us to solve the motion of the host fluid with periodic 
boundary conditions. Couette flow profiles are approxi- 
mately formed over a range from y = —L y /A to y = L y /A, 
although the zigzag profile strictly becomes a Couette 
flow when L y — > 00. <j)f p represents the body force that 
ensures the rigidity of particles and the aprropriate non- 
slip boundary conditions at the fluid/particle interface, 
which is further elaborated in reference [l(| [0, EH • 

III. RESULTS AND DISCUSSION 

The computational domain has three dimensions 64 x 
64 x 64 and periodic boundary conditions. The numeri- 
cal parameters for both the host fluid and the spherical 
particles are A = 1, 77 = 1, pf = 1, a = 5, £ = 2, 
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and particle density p p = 1. The imposed shear rate is 
7 = 0.005, and the temperature is fc^T = 0.07. The tem- 
perature was determined by measuring the long-time dif- 
fusion coefficient of a single Brownian particle at thermal 
equilibrium before the shear was imposed. This system 
has dimensionless parameters such that the Peclet num- 
ber Pe = 6iri]a 3 j /ksT ~ 170 and the particle Reynolds 
number Re p = p p A fa 2 /ii = 0.125. The initial configu- 
ration of the spherical particle is located at the central 
position of the system. The volume fraction of a single 
particle is $ = 0.002. 

An important feature of a single Brownian particle in 
a simple shear flow is that the MSD in the flow direc- 
tion varies asymptotically with time as t 3 . Although we 
would expect t-dependence of the MSD in the long-time 
limit, this non-diffusional behavior is explained as a result 
of a coupling between diffusive motion in the velocity- 
gradient direction and convective motion due to shear in 
the flow direction. Theoretical solutions to this problem 
show the MSD in the flow direction^) is composed of 
three parts: 

(RUtf) = vl(yo)t 2 + (RUt) 2 h + (Ri(tf)o (6) 

assuming that, at t = 0, the particle is in the initial 
position (0, yo, 0), wherein the first term on the right 
hand side of Eq. (J6J) is a shear contribution only, and 
represents the simple translation of a particle along the 
streamline in the flow direction. The second term is the 
coupling term described above, which represents the con- 
vection induced by diffusion. The third term, which is in 
the same form as the MSD of a free Brownian particle, 
is a thermal contribution only. Similar to that of the x 
direction, the MSDs in the velocity-gradient direction(y) 
and the vorticity direction(z) are derived theoretically; 
{R\{tf) = (i?f(i) 2 ) and <i?f(i) 2 } = (Rtm . These 
quantities arc determined analytically for a single Brow- 
nian particle obeying the normal Langevin equation (LE) 
or the generalized Langevin equation (GLE) with mem- 
ory effects [IE 03, El- 

The time evolution of the MSD in the x direction (Q>) 
for a single spherical particle under shear is shown in 
Fig. [2 The MSD was calculated via (\Rf(t) - Rf (0) - 
v x (yo)t\ 2 ) to eliminate the purely shear contribution that 
corresponds to the first term of Eq. ©. The MSD is 
scaled by ksT. The MSD is in excellent agreement with 
the analytical solution for the GLE. For short times of up 
to t ~ 10 2 , the motion of the Brownian particle in the x 
direction is like that of a free Brownian particle obeying 
the GLE. For long times, t > t„ = 25, the MSD asymp- 
totically approaches 2Doj 2 t 3 /3, where Do = ksT /Qivqa. 
This t 3 regime in the MSD is approached in a much slower 
manner, t~ 1 ' 2 , than that of the analytical solution for the 
LE. The transient behavior of the MSD depends strongly 
on whether the memory effects are taken into account. 

The MSD (\Rf(t) - Rf(0)| 2 ) in the z direction (■) 
is also plotted in Fig. [5] The numerical results of the 
MSD in the z direction agree well with the analytical 
solution for the GLE of a free Brownian particle where 



the t regime is approached in a much slower manner, 
i -1 / 2 , and the diffusive motion is attained on time scales 
of O(10 4 ). The diffusion time characterizing the diffusive 
motion is td = a 2 I Do ~ 3.4 x 10 4 , which measures the 
particle diffusion over the particle size. 

Another dynamical quantity of interest is the posi- 
tional cross-correlation (i?f (t)R^ (t)) of a Brownian parti- 
cle in a simple shear flow. This analytical form is derived 
for both the LE and GLE. These long-time behaviors 
show the same dependence on t, which is Do^t 2 . The 
cross-correlation was calculated via ((Rf(t) — Rf{0) — 
Vf{0)t){R?{t) - Rf(0))). The numerical results (A) arc 
plotted in Fig. [5] and are in good agreement with the 
analytical solution for the GLE. 

This detailed analysis of the MSD is sufficient to con- 
firm the validity of our method of incorporating memory 
effects; however, within our method, it is assumed that 
the thermal fluctuations can be represented by white and 
gaussian noise. This means that correlations between 
thermal noise at different times are completely ignored, 
although in the theoretical framework of the GLE, noises 
are correlated. The thermal noise memory, which we ig- 
nore, is considered to affect the particle's motion on very 
fast time scales t < tb ~ 4.5. On these very short time 
scales, t < tb, a gap between the hydrodynamic analyti- 
cal solutions and the numerical results is clearly observed 
in previous studies [HQ3. Although our method is not 
applicable to particle motion for t <tb, the method can 
be accurately applied to the motion of particles in a shear 
flow on time scales t > tb, in which the memory effects 
become significant. 

Most numerical approaches so far performed for con- 
centrated particles dispersions are based on the Stokcsian 
dynamics thus those simulations are valid only for 
the zero Reynolds number limit. Recently, concentrated 
dispersions at finite Reynolds number have been simu- 
lated using the DNS methods based on the stochastic 
rotation dynamics po| or the lattice Boltzmann method 
[HI, HI We have applied the modified SPM to 

concentrated dispersions at a finite Reynolds number 
Re p = 0.125 in order to examine the short-time motion 
of Brownian particles in a shear flow at finite volume 
fractions. Simulations were performed at several volume 
fractions from a very dilute case with $ = 0.002 to very 
dense cases with $ = 0.4 and 0.5, which arc apparently 
higher than the previous DNS simulations. 

Since Pe ~ 170, the hydrodynamic shear forces be- 
come dominant over the thermal forces in the particle' 
motion. The initial configurations of the dispersed par- 
ticles are randomly distributed. 

Figure [3] shows the time evolution of the MSD in the 
z direction for several volume fractions, up to $ = 0.5. 
The time is scaled by I/7. As the volume fraction is in- 
creased, the MSD grows more rapidly in time due to the 
hydrodynamic and direct interactions between the parti- 
cles, resulting in a increase in the slope of the MSD for 
a small strain jt ~ O(10 _1 ). For high volume fractions 
0.2 < $ < 0.4, the diffusive behavior is attained at a 
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FIG. 2: Mean square displacement scaled by teT of a sin- 
gle spherical particle fluctuating in a Newtonian fluid under 
shear at k B T = 0.07 and 7 = 0.005: (O) now direction, (■) 
vorticity direction, and (A) positional cross-correlation. The 
analytical solution of the MSD in the flow direction (x) for a 
Brownian particle in a shear flow: the dashed line represents 
the LE and the solid line the GLE. The dash-dotted line is 
the analytical solution of the MSD in the vorticity direction 
(z) for the GLE, which is the same form as that for the GLE 
of a free Brownian particle The analytical solution of the po- 
sitional cross-correlation for a Brownian particle in a shear 
flow: the dotted line represents the LE and the dashed line 
the GLE. The Brownian time tb — 4.5, the kinematic time 
tv = 25, and the diffusion time to — 3.4 x 10 4 . The parame- 
ters for both the particle and fluid are a = 5, £ = 2, p p = 1, 
pf — 1, and 77 = 1. The time unit is A 2 jv = 1, and the length 
unit is A = 1. 



smaller strain yt ~ 0(1) rather than at <t> = 0.002, as for 
a Brownian particle. The accelerated Stokcsian dynamics 
(ASD) for non-Brownian particles at $ = 0.2 show the 
diffusive region is attained at larger strains than at least 
10 [HJ Compared with these numerical results, the onset 
of the diffusive region in the present simulations is much 
faster. The time at which the diffusive region is attained 
shifts to shorter times at higher volume fractions. 

For the highest volume fraction <t> = 0.5, however, we 
see that the motion of Brownian particles are trapped 
within the effective cages formed by the surrounding par- 
ticles. This is because the particles start to form string- 
like objects in the flow direction, and finally the whole 
system evolves into a two-dimensional ordered structure 
in the plane perpendicular to the flow, similarly to the 
flow-induced ordering commonly observed in experiments 
[25l l2rl l27t and by computer simulations [H, |29, [3(| un- 
der shear flow. In the present case with Re p = 0.125 
and Pe ~ 170, we found that the shear-induced ordering 
occurs at the high volume fraction between <!> = 0.4 and 
0.5. 



The long-time diffusion coefficient D z was calculated 
by linearly fitting the data over the diffusion regions. The 
inset of Fig. [3] shows the volume fraction dependence 
of D z scaled by 7a 2 . The diffusion coefficient increases 
rapidly up to <!> = 0.3 and reaches a plateau at a volume 
fraction beyond $ > 0.3. This behavior of D z is remark- 
ably different from that at thermal equilibrium, where 
the diffusion coefficient decreases with increasing volume 
fraction. The enhancement of the diffusion coefficient 
with increasing volume fraction is a typical characteris- 
tic of shear-induced diffusion coefficients. These results 
exhibit the same qualitative behavior as the experimental 
results obtained for non-Brownian particles [3l|, HH, [33| . 
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FIG. 3: The time evolution of the mean square displacement 
in the vorticity direction (z) for several volume fractions $ at 
k B T = 0.07 and 7 = 0.005. The Peclet number Pe ~ 170 and 
the particle Reynolds number Re p = 0.125. The inset shows 
the volume fraction dependence of D z scaled by 7a 2 . The 
parameters for both the particle and fluid are a — 5, £ = 2, 
p p — I, pf = 1, and 77 = 1. 



IV. CONCLUSION 

In conclusion, the short-time motion of Brownian par- 
ticles in an incompressible Newtonian fluid under shear 
was examined by using a modified SPM in which exter- 
nal forces are introduced to approximately form a simple 
shear flow throughout the entire system with periodic 
boundary conditions. The validity of the method was 
carefully examined by comparing the present numerical 
results for the MSD with the hydrodynamic analytical 
solution of the GLE of a single Brownian particle in a 
simple shear flow. In the present study, we aimed to 
modify the original SPM by incorporating thermal fluc- 
tuations so that the modified SPM is valid for tb < t, 
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while other computational methods such as Brownian 
dynamics and ASD, which are based on the LE that 
ignores memory effects due to fluid motions, are valid 
only for rg <C td < t. Simulations were then performed 
for monodisperse dispersions of repulsive spherical par- 
ticles at volume fractions ranging from $ = 0.002-0.50. 
We found that the MSD in the vorticity direction grows 
rapidly in time and with increasing particle volume frac- 
tion. At a strain jt ~ 0(1), the diffusive region is at- 



tained for 0.2 < $ < 0.4. The onset of the diffusive 
region shifts to shorter times at higher volume fractions. 
For $ = 0.5, however, the particles are no more diffu- 
sive because of the shear-induced ordering. The diffusion 
coefficient in the vorticity direction was obtained from 
the long-time behavior of the MSD. For volume fractions 
of up to 0.3, the diffusion coefficient rises rapidly with 
increasing volume fraction. It then levels off for volume 
fractions beyond 0.3. 
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